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Abstract 

Background: Planl<tonic ciliated larvae are characteristic for the life cycle of marine invertebrates. Their most 
prominent feature is the apical organ harboring sensory cells and neurons of largely undetermined function. An 
elucidation of the relationships between various forms of primary larvae and apical organs is key to understanding 
the evolution of animal life cycles. These relationships have remained enigmatic due to the scarcity of comparative 
molecular data. 

Results: To compare apical organs and larval body patterning, we have studied regionalization of the episphere, 
the upper hemisphere of the trochophore larva of the marine annelid Platynereis dumerilii. We examined the spatial 
distribution of transcription factors and of Wnt signaling components previously implicated in anterior neural 
development. Pharmacological activation of Wnt signaling with Gsk3|3 antagonists abolishes expression of apical 
markers, consistent with a repressive role of Wnt signaling in the specification of apical tissue. We refer to this 
Wnt-sensitive, six3- and foxc^2-expressing part of the episphere as the 'apical plate'. We also unraveled a molecular 
signature of the apical organ - devoid of six3 but expressing foxj, irx, nkx3 and box - that is shared with other marine 
phyla including cnidarians. Finally, we characterized the cell types that form part of the apical organ by profiling by 
image registration, which allows parallel expression profiling of multiple cells. Besides the /lox-expressing apical tuft 
cells, this revealed the presence of putative light- and mechanosensory as well as multiple peptidergic cell types 
that we compared to apical organ cell types of other animal phyla. 

Conclusions: The similar formation of a six3+, foxq2+ apical plate, sensitive to Wnt activity and with an apical tuft 
in its 5/x3-free center, is most parsimoniously explained by evolutionary conservation. We propose that a simple 
apical organ - comprising an apical tuft and a basal plexus innervated by sensory-neurosecretory apical plate 
cells - was present in the last common ancestors of cnidarians and bilaterians. One of its ancient functions would 
have been the control of metamorphosis. Various types of apical plate cells would then have subsequently been 
added to the apical organ in the divergent bilaterian lineages. Our findings support an ancient and common origin 
of primary ciliated larvae. 
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Background 

The vast majority of animal phyla live in the ocean and 
develop via small ciliated larvae that form part of the 
zooplankton [1]. These larvae, called primary larvae', are 
equipped with sensory cells to perceive various stimuli in- 
cluding light, touch and chemical cues [2,3]. Simple ner- 
vous systems integrate sensory information and control 
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ciliary locomotion [4]. Most conspicuously, an apical 
organ' is found in various groups, including cnidarians 
[5], protostome annelids, mollusks [2], flatworms [6] and 
nemertines [7], as well as deuterostome echinoderms [8] 
and hemichordates [9]. These groups belong to the 
Neuralia (which includes cnidarians and bilaterian proto- 
stomes and deuterostomes; Figure 1) [10]. Apical organs 
are often equipped with an apical tuft' of long cilia and in- 
clude a small set of sensory-neurosecretory cell types 
[11,12]. As shown for cnidarians, mollusks and annelids, 
apical organs play a role in the control of settlement 
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Cnidarian Lophotrochozoan Deuterostome 



pelagic benthic pelagic benthic pelagic benthic 




/\ Pelago-benthic neuralian ancestors: Single terminal addition 
of benthic adult stages (Jagersten 1972; Arendt et. al, 2013) 




B Holopelagic neuralian ancestors.: Terminal addition of benthic 
adult benthic stages two (or three) times (Nielsen 2012) 




Q Holobenthic neuralian ancestors: Repeated intercalation of 
swimming larval stages; many times in each lineage (Raff 2008) 



f ^ 

Figure 1 The three scenarios show the origin of the pelagic 
larval body plan, indicated by red arrows. The presence of 
pelagic forms is indicated by blue lines and that of benthic forms by 
brown lines. A single-coloured line indicates a monophasic life 
cycle that would be pelagic in scenario B and benthic in scenario C. 
Double lines (blue and brown) indicate a biphasic, pelago-benthic 
life cycle (with pelagic larval and benthic adult forms). Note that the 
biphasic life cycle is assumed to have evolved multiple times 
independently in scenario C. (A) The classical view implies 
homology of both ciliated larvae and benthic adults that, once 
evolved, have remained part of the eumetazoan life cycle [17]. 
(B) Nielsen [16] modified this view to propose that the holopelagic 
neuralian ancestors persisted beyond the initial divergence of the 
major neuralian clades, and that the biphasic life cycle with benthic 
adults arose independently in the cnidarians and once or twice in 
the bilaterians. (C) In stark contrast, other authors assume that 
today's ciliated larvae arose convergently many times by the 
repeated intercalation of a pelagic dispersal larva into primarily 
monophasic, holobenthic life cycles and are thus evolutionarily 
unrelated. This view implies that the characteristics of today's 
swimming larvae such as apical organs and equatorial ciliary bands 

evolved convergently [18]. 

I. J 

[13-15], facilitating the integration of multimodal sensory 
input and the coordination of effector cells. Development 
involving primary larvae (or any other kind of secondary 
larvae) is referred to as Indirect'; if it does not involve a 
larval stage, it is referred to as 'direct' [16]. 

Widespread occurrence of ciliated larvae, with similar 
overall body plans, has prompted radically different views 
regarding their evolutionary significance [16,18-20]. While 
some authors consider them an ancient feature of the 
metazoan life cycle [16,17,19], others assume that they 
evolved multiple times independently [18,20] (Figure 1). 
In the past two decades, studies of developmental genes 
involved in setting up the larval body plan have started to 
provide a new source of comparative data to resolve these 
conflicting views. These studies have revealed specific 
similarities such as the expression of the goosecoid and 
brachyury genes in foregut, midline and hindgut of lopho- 
trochozoan and deuterostome larvae [21] (in a manner 
distinct from direct developers such as insects or verte- 
brates). Also, the ciliary bands characteristic for swimming 
larvae were found to specifically express otx and three 
conserved microRNAs in lophotrochozoans and deutero- 
stomes, corroborating the homology of shared larval fea- 
tures [22,23]. Conversely, dissimilar expression reported 
for nk2,l and hnfS in the apical plate of sea urchins and 
the episphere of abalone larvae [24] raised doubt regard- 
ing the common origin of similar shared larval features. 
However, these pioneering studies thus far relied on small 
gene sets. 

More recently, an extended set of transcription factors 
has been shown to play a conserved role in patterning 
the larval body plan in hemichordate [25] and in 
sea urchin [26,27], belonging to the deuterostomes, and 
in the sea anemone Nematostella vectensis, a cnidarian 
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[28]. These factors respond to the differential stabilization 
of |3-catenin along the primary body axis, triggered by 
Wnt signaling [28,29]. In particular, six3 and foxq2 have 
been shown to negatively respond to Wnt signaling in a 
complex sequence of patterning events [28,29] . These fac- 
tors specify apical territory around the apical pole [27,30] . 
For the first time, the conserved regional expression of 
similar sets of transcription factors provides a molecular 
framework for the comparison of larval cell types and 
tissues and thus important clues to larval body plan and 
apical organ evolution. 

In our current study, we investigated the apical pattern- 
ing system in the marine annelid Platynereis dumerilii 
[31], a lophotrochozoan protostome with a canonical bi- 
phasic pelago-benthic life cycle (involving a pelagic phase, 
that is, swimming, primary larvae, and a benthic phase, 
that is, bottom-dwelling adults). We examined the set of 
transcription factors involved in apical patterning, includ- 
ing six3 and foxq2, in the episphere of the Platynereis 
trochophore larva and have shown that, as in deutero- 
stomes and cnidarians, expression of these factors is sensi- 
tive to Wnt signaling. We found that the apical organ 
develops in a small central territory devoid of sixS expres- 
sion, that instead expresses a number of other factors, 
many of which are found in the same location in other 
neuralians. By expression profiling, we molecularly char- 
acterized several cell types that form part of the apical 
organ in Platynereis, which we compared to apical organ 
cell types described for other animal groups. 

Our results reveal that the larvae of cnidarians, proto- 
stomes and deuterostomes exhibit extensive similarity in 
the molecular topography of body regions around the ap- 
ical organ, which we use to genetically define apical plate' 
and apical organ'; the specification of these regions by a 
conserved apical signaling system; and the molecular fin- 
gerprint of a subset of apical organ cell types. These find- 
ings support homology of some primitive type of apical 
organ (and thus of swimming ciliated larvae) in Neuralia 
and are most consistent with an early and unique origin of 
animal larval forms (Figure 1A,B). 

Results 

Molecular topography of the apical region 

In the sea urchin [26,27], the hemichordate Saccoglossus 
kowalevskii [25] and the anthozoan cnidarian Nematostella 
[28,30], expression of the transcription factors sixS and 
foxq2 demarcates the most apical body region. In Saccoglos- 
sus and Nematostella, the six3-\-,foxq2-\- domain is peripher- 
ally overlapping with a ring of rx and bounded by even 
more peripheral rings of otx and otp expression [25,28]. 
Likewise, is expressed in the apical plate in sea urchin 
[32] and the hemichordate [33]. We have previously 
shown that six3 is expressed in a large contiguous domain 



of the Platynereis episphere [34], peripherally overlapping 
with the expression of rx [35]. 

Building on this, we set out to further refine the mo- 
lecular topography of the Platynereis larval episphere. 
We found that, in Platynereis, the newly characterized 
foxq2 expression demarcated the upper' two thirds of 
the episphere (Figure 2A), largely overlapping with that 
of six3 (Figure 2B). We refer to this six3+, foxq2-\- region 
as apical plate, in accordance with classical nomenclature 
for spiralian trochophore larvae (for example, see [2] and 
discussion). Expression of the newly characterized fezf 
(Figure 2C) and that of rx (Figure 2D) overlapped that of 
six3 and foxq2 peripherally (with rx being largely restricted 
to the dorsal and fezf to the ventral body sides). More per- 
ipherally, in the lower' third of the episphere, we detected 
expression of otx (Figure 2E) [21], partially overlapping 
with that of rx dorsally and that of fezf ventrally. Even 
more peripherally, expression of the otp gene (Figure 2G) 
[12] demarcated the ciliary band, called the prototroch. 

The Platynereis episphere can thus be subdivided into 
a sequence of molecular regions, arranged in concentric 
rings from apical to peripheral (Figure 3). This overall 
sequence matches the molecular topography observed in 
deuterostomes and cnidarians (Figure 3). Another con- 
served feature of apical patterning in sea urchin [24,36], 
Saccoglossus [25] and again in Platynereis [12] is expres- 
sion of nkZl in an apical-ventral territory (Figure 2F), 
overlapping six3. In sea urchin, this domain includes 
cells positive for onecut {HNF6) [24], and we found that, 
in Platynereis, expression of the same genes covers cells 
located ventromedially at a similar distance to the apical 
tuft cells (demarcated by a yellow dotted line in the con- 
focal reflection microscopy panels; compare Figure 2F" 
with Figure 2J'), which indicates co-expression. This is im- 
portant because the two genes were shown to comprise 
non-overlapping domains in the apical region of the red 
abalone Haliotis (where onecut is expressed dorsal to the 
apical organ only [24]; note that Platynereis also shows 
dorsal patches of onecut expression in Figure 2J'). 

The Platynereis trochophore also shares expression 
features with cnidarian and deuterostome larvae spe- 
cific to the developing apical organ. In sea urchin [37], 
Nematostella and Platynereis, the apical tuft and a girdle 
of cells surrounding the tuft are devoid of six3 expression 
(white dashed line in Figure 2B' and compare to Figure 3 
in Poustka et al. [37] and Figure 3d in Marlow et al [28]). 
In Platynereis, the six3-ivQe: region extends into a dorsal 
patch of cells that becomes separated from the central 
patch by a band of six3^- cells (which includes the con- 
spicuous crescent cells, see below). Inside the six3-ivQQ re- 
gion, we detect restricted expression of several 
transcription factors (using the nuclei of the two apical 
tuft cells for reference; yellow circles in confocal apical 
views in Figure 2): otp (Figure 2G"; in addition to its ring- 
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(See figure on previous page.) 

Figure 2 Apical regional patterning mRNA expression in early trochophore larvae. (A-G) are ventral views, (A'-G') and (J-Q) are apical 
views, and (A"-G") and (J'-Q') are apical views using confocal reflection microscopy of NBT/BCIP staining (red) following in situ mRNA and 
fluorescent DAPI staining of nuclei (white), allowing localization of individual cells in the apical Plate. F", I, J', K', L', P' and Q' are high magnification 
apical views. Ampullary cells are demarcated with yellow dashed circle. All images are of 24 hour post-fertilization (hpf) embryos, except C, C, D, 
D', E, E', G, G', L, M, 0, 0', which are of 20 hpf embryos. Asterisks mark the tip of the apical pole. (A) Broad apical expression of foxq2 in 24 hpf trocho- 
phores. (B) Broad apical expression of s/x3, excluding area encompassing apical organ and a dorsal patch of cells (demarcated with white dashed lines). 
(C) Fezf ventrally in apical plate. (D) Rx in dorsal and ventral domains between the six3 domain and the prototroch. (E) Otx in 
ring-like domain adjacent to prototroch. (F) A//c2./ in a ventral strip of cells and in two spots flanking the crescent cells. (G, I) Otp in apical cluster and 
in prototroch. (H, H') Schematized ventral and apical views at 24 hpf depicting crescent and ampullary cells and prototroch. (J) Onecut in crescent cells 
and in cells ventral to the ampullary cells. (K) FoxJ in the ampullary cells, crescent cells and prototroch. (K') Antibody staining to acetylated tubulin 
(green) and DAPI-labeled nuclei (blue). (L) Hoxl in ampullary cells. (M) Irx in s/xi-negative territory and in an equatorial ring below the prototroch. 
(N) FgfR in a cell dorsal to the ampullary cells. (O) Noggin in a cell dorsal to the ampullary cells. (P) trpV\n cells ventral to the ampullary cells. (Q) Nl<x3 
in a small patch of cells ventral to the ampullary cells. (Q') Nl<x3 in mechanoreceptor cells (yellow arrows). 



like expression, see above), nkx3 (Figure 2Q'), foxj (Fig- 
ure 2K') and hoxl (Figure 2U). In addition, the irx gene is 
expressed in the dorsal patch of cells devoid of six3 
expression (Figure 2M'). Notably, nkx3 [28], foxJ [30], the 
hox gene paralog anthoxl [28] and irx [28] are also 
expressed in the apical region devoid of six3 in Nematos- 
tella; foxj matches the region devoid of six3 expression in 
sea urchin [38] and apical irx expression has been docu- 
mented in hemichordates, although broader than in anne- 
lids [25] . An apical region devoid of six3 expression. Tilled' 
by the restricted expression of foxj, nkx3, hox and irx, thus 
represents a larval-specific transcription factor signature 
that molecularly defines the apical organ in several dis- 
tinct, phylogenetically remote neuralian groups such as 
anemone and annelid (Figure 3) and lends strong support 
to its evolutionary conservation (see Discussion). 

Functional studies have demonstrated the role of fibro- 
blast growth factor (Fgf), transforming growth factor 



(Tgf|3 and Fgf signaling) and Wnt signaling in apical plate 
patterning. For example, Fgf signaling controls formation 
of the apical organ in sea anemone [14]; the apical plate of 
sea urchins is patterned by Tgfp signaling [39] and fgfrl 
localizes to the apical plate during apical organ formation 
[37]. Similarly, the Platynereis apical organ region showed 
specific expression of fgfr in a highly restricted population 
of cells (Figure 2N). We also found an apically restricted 
domain of the Tgfp signaling antagonist noggin in cells de- 
void of six3 expression (white dashed line in Figure 20'), 
which indicates that both signaling systems also play a 
role in apical patterning in Platynereis, Notably, ex- 
pression of noggin has also been reported for the six3' 
free apical organ in sea anemone [40]. Consistent with a 
prominent role of Wnt signaling in Platynereis apical pat- 
terning, we detected the Wnt receptor mo\QC\AQ frizzledS/ 
8 in a broad apical domain (Figure 4A; as seen in the 
hemichordate Ptychodera [41] and sea urchin [42]) and 




Figure 3 Comparison of apical molecular territories in a sea anemone planula larva and an annelid trochophore. Apical (above) and lateral 
(below) views of a schematized planula larva of a sea anemone (left panel). Gene expression based on published data for Nematosteiia vectensis (for 
references see text). Apical (above) and lateral (below) views of a schematized annelid trochophore larva (right panel). Gene expression based on our data. 

V J 
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Figure 4 Wnt pathway members are expressed in the episphere and hyposphere during early development. The prototroch is marked 
with a dashed line. The tip of the apical plate is marked with an asterisk. Arrows indicate dorsal (D), ventral (V), anterior (A) and posterior (P). 
(A) Frizzled5/8 across the apical plate, 16 hours post-fertilization (hpf) apical view and a (A') 20 hpf ventral view (inset shows a surface ventral 
view). (B) Wnt4 along the blastopore, 16 hpf lateral view. (B') Confocal image of ventral view of Wnt4 expression in which antibody staining to 
acetylated tubulin is visible with a green secondary antibody, reflection signal of /n situ hybridization is red and nuclei are visible, stained with 
DAPI, in blue. (C) Apical view of sfrp 1/5 expression at 24 hpf and at (C) 48 hpf. (A", B" and C") Schematic representations of the expression pat- 
terns in (A) and (B) and (C) and (A') and (B') and (C). 



sfrp, which plays a conserved role in antagonizing Wnt 
signaling in deuterostomes [42,43], specifically in the ap- 
ical organ region (Figure 4C,C',C")- Furthermore, we de- 
tected early expression of wnt4 in ventral peripheral 
episphere regions, outside of (but abutting) the six3 
domain, consistent with a role in apical patterning 
(Figure 4B,B'). Wnt4 is one example of a Wnt paralog 
expressed at the right stage and location in developing 
Nematostella [44] and amphioxus [45] to be involved 
in patterning along the primary axis; yet, in the ab- 
sence of specific knockdown or knockout data, no case 
for the specific involvement of Wnt4 can be made. 
Note that none of the many Wnts expressed in Platynereis 
larvae [46,47] appears to be expressed in the sixS^- apical 
plate, in line with repression of Wnt signaling being a pre- 
requisite for apical specification to occur, a notion that we 
set out to test further experimentally. 



Ectopic Wnt signaling abolishes specification of the apical 
organ region 

In deuterostome larvae such as sea urchin and Saccoglossus, 
as well as in the cnidarian Nematostella, the expression of 
apical markers such as foxq2 and sixS critically depends on 
levels of active Wnt signaling [27,43,48]. Activation of ca- 
nonical Wnt signaling via chemical inhibition of Gsk3|3 
via Li + and/or paullone treatment in sea urchins [37], 
hemichordates [43] and cnidarians [28] results in the loss 
of apical ectodermal markers. Conversely, degradation of 
nuclear p-catenin in sea urchin leads to expansion of 
apical ectodermal genes including nk2,l, foxq2, sixS, rx 
and fgfR [27], indicating that Wnt- mediated antagonism 
of apical plate markers plays a role in the development of 
apical territories in deuterostome and cnidarian larvae. 
To test the role of Wnt signaling in apical patterning 
in Platynereis larvae, we exposed early trochophores to 



Marlow et al. BMC Biology 2014, 12:7 
http://www.bionnedcentral.conn/1 741 -7007/1 2/7 



Page 7 of 1 7 



foxQ2 



hox1 



CD 





Figure 5 Azakenpaullone treatment of embryos alters gene expression in the episphere in a dose-dependent manner. All images are 
apical views with dorsal at the top of the image and ventral at the bottom. (A) Six3 is severely reduced in embryos treated with azakenpaullone. 
(B) Foxq2 expression is lost in all but the most apical cells of the episphere in azakenpaullone treatments. (C) Fewer otp cells are stained at 
increasing concentrations of azakenpaullone and expression is entirely lost at 10 pM. (D) Hoxl, a molecular marker of the ampullary cells (arrows), 
is not affected by azakenpaullone treatment. Ampullary cells are specified very early in development, earlier than the initiation of chemical 
treatments at 12 hours post-fertilization, prior to the differentiation of the remainder of the episphere. (E) Pax6 expression, which marks the 
ventral regions of the developing apical plate, is expanded in the episphere upon treatment with azakenpaullone. Counts of the embryos 
displaying wild-type, reduced or expanded expression domains for two replicates at 0, 1 and 5 concentrations are displayed in Figure S4. 



azakenpaullone, a selective inhibitor of GskSp [49] that 
has been shown to trigger nuclear p-catenin accumulation 
in Platynereis, mimicking ectopic activation of Wnt sig- 
naling [47,50]. We found that azakenpaullone specifically 
knocked down or abolished apical expression of episphere 
markers, including the broadly expressed sixS (Figure 5A) 
and foxq2 (Figure 5B) in a concentration-dependent man- 
ner, with the majority of expression reduced or lost be- 
tween 1 (iM and 5 (iM (Figure 5). Conversely, the 
expression of pax6, which occupies a more ventral periph- 
eral position in the larval episphere [51] (overlapping with 



Wnt4; Figure 4E), was expanded at 0.5 (iM to 10 (iM 
concentrations of azakenpaullone (Figure 5E). The num- 
ber of otp + apical organ cells was also reduced at higher 
concentrations of azakenpaullone (Figure 5C), with all 
cells absent at 5 and 10 (iM concentrations. By contrast, 
the expression of hoxl in the apical tuft cells persisted 
at all concentrations (Figure 5D, arrows). This may 
have been because tuft cells are among the first cells 
to differentiate in the apical plate (Figure 6A). To 
test the dynamic role of Wnt signaling in episphere 
patterning, we conducted washout experiments in 
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Figure 6 Molecular fingerprint of Platynereis apical organ cell types. (A) EdU incorporation (white) in dividing cells between 22 and 
24 hours post-fertilization (hpf) (apical view). Central EdU-negative cells (magenta) represent post-mitotic apical organ cells. Inset: 5 |jm z- 
projection of apical organ. (B) Acetylated tubulin staining at 30 hpf, apical view. Crescent cells located dorsal to ampullary tuft cells. Inset: Cilia of 
the putative mechanosensory cells (arrowhead). (C) Lateral view of ampullary tuft cells, acetylated tubulin (green) and DAPI staining (blue), 30 hpf 
(D) Serotoninergic cells of the apical organ stained with an anti-5HT antibody, 48 hpf. White arrow indicates serotonergic interneuron (magenta 
in schematic). (E) Four otp + peptidergic flask-shaped cell (numbers) at 48 hpf (apical view). One of these labeled by an FMRFamide antibody. In 
schematic, numbers indicate otp + cells and the FMRFamidergic cell is in cyan. (F) Lateral view of putative mechanosensory cells, 48 hpf (arrow- 
heads: stiff curly tips of the sensory cilia). (G, H) Mutually exclusive gene expression at 48 hpf, 23 pm z-projection of averaged expression patterns 
obtained with profiling by image registration (PrImR). (I) Cell types of the otp + domain. The image shows non-overlapping expression of markers 
of each cell as deduced from PrImR, 1 1 pm z-projection of average expression patterns. (J) 4 pm z-projection of the average axonal scaffold at 48 
hpf, showing the deep position of the two ampullary cells surrounded by neuropil. (K) Expression of peropsin at 24 hpf (L) Expression of c-opsinl 
at 20 hpf (M) Hierarchical clustering of the molecular fingerprint of individual cell types identified within the apical organ region and of the pro- 
totroch for comparison. Expression in magenta. Asterisk indicates the two ampullary tuft cells, red arrows indicate crescent cells and yellow arrow- 
heads indicate cilia of two putative mechanosensory cells. In panels B-H, gene expression is in red, DAPI staining in blue and tubulin staining in 
green unless otherwise specified. Panels C-F, scale bar is 10 pm. 



azakenpauUone-treated embryos (Additional file 1: 
Figure S3). Following washout at 24 hours, embryos 
were assessed at 30 hours post-fertilization (hpf) for 
recovery of gene expression. We saw moderate recov- 
ery in the expression of foxq2 and sixS as well as a 
slight restriction of pax6, but no change in the number 
and location of otp -h cells. We attribute the moderate re- 
covery in expression to the determinate lineage of Platy- 
nereis larvae, in which stereotyped divisions may result in 
a restricted fate potential very early in development. Taken 
together, our data indicate that in Platynereis as in deu- 
terostome larvae, the transcription factors defining the 
molecular identity of apical body regions, of which sixS 



and foxq2 form core components, are opposed by a Wnt- 
dependent signaling center. 

Morphological and molecular characterization of apical 
organ cell types 

Light and electron optic studies have revealed cell types 
that make up the apical organs of different invertebrate 
larval groups [11,52-54]. In Platynereis^ the apical tuft 
was visible by 16 hpf (cf. ^r/?F channel expression at 24 
hpf. Figure 2P) and other apical organ cells were likewise 
post-mitotic by 24 hpf (Figure 6A), expressing markers 
indicative of neuronal differentiation. At 30 hpf, several 
apical cell types were fully morphologically differentiated 
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(Figure 6B; [31]). We found that the apical tuft was 
formed by two basket-shaped cells with intracellular 
tubulin support structures (Figure 6C); cells with a very 
similar morphology, referred to as ampullary cells, have 
previously been described in moUusk larvae [11]. These 
cells persisted deep in the medial brain at later stages in 
the center of a massive commissural and neurosecretory 
neuropil (Figure 6J), and may thus represent a structural 
organizing center for the juvenile nervous system, as 
suggested for other polychaete larvae [50]. Dorsal to the 
ampullary tuft cells, we found another set of large cells 
with multiple motile cilia in a crescent-moon shape, 
known as crescent cells (Figure 6D). Two serotonergic 
cells have also been found in the apical organ region 
by 30 hpf [31]. Closest to the tuft was a serotonergic inter- 
neuron (white arrow. Figure 6D) lacking sensory den- 
drites. This cell was located deep in the epithelium, 
adjacent to an assembly of previously described sensory- 
neurosecretory flask-shaped cells (Figure 6E) that, mor- 
phologically, resemble chemosensory cells [12] (called 
parampullary cells in moUusks [11]). More ventral to the 
parampullary cells, we detected a median pair of cells 
bearing short, stiff and curly sensory cilia resembling 
mechanoreceptors (Figure 6F). 

These distinctive morphologies, in conjunction with 
the recently established Profiling by Image Registration 
(PrImR) technique, enabled us to assign a molecular fin- 
gerprint to these cells, providing them with unique mo- 
lecular identities [47]. PrImR utilizes the stereotyped 
development of the Platynereis axonal scaffold to gener- 
ate in silico alignments of mRNA in situ expression pat- 
terns, and allows single cell co-expression analyses to be 
conducted [47]. Of a collection of 140 genes currently 
available for PrImR single cell co-expression analysis, 29 
were differentially expressed in cells of the apical organ 
region in the 48 hpf larva. Additional file 1: Figure S5 
details the PrImR-based co-expression analysis for the 
above-mentioned morphologically identifiable cells, 
namely the ventral-most serotoninergic cell (Additional 
file 1: Figure S5A), the parampullary sensory-neurosecretory 
cells (Additional file 1: Figure S5B-E), the ampullary tuft 
cell (Additional file 1: Figure S4F), the crescent cells 
(Additional file 1: Figure S5G) and the pair of puta- 
tive mechanoreceptors (Additional file 1: Figure S5H). 

PrImR revealed unique sets of genes expressed by each 
of these cell types, in line with their specialized sensory- 
neurosecretory and neuronal characteristics. For example, 
as determined previously, the flask-shaped parampullary 
cells expresses otp (Figure 6G; Additional file 1: Figure 
S5B-E; and see above), mir-7 (Figure 6H) and prohormone 
convertase 2 {phc2) [12]. Beyond that, PrImR allowed 
cellular allocation of transcripts encoding neuropep- 
tide precursors for DLamide, FMRFamide and WLDa- 
mide (Figure 61), consistent with a conserved role of 



otp in specifying different types of peptidergic cells 
[12]. These three neuropeptides mediate opposing effects 
on locomotor behavior with DLamide and FMRFamide 
exposure leading to an increase and WLDamide leading 
to a decrease in ciliary beating frequency [55]. The otp + 
peptidergic cells also expressed MIP, the recently described 
settlement-inducing neuropeptide [15]. In addition, we 
found that all otp + sensory-neurosecretory cells were posi- 
tive for the newly identified peropsin gene (Figure 6K), an 
opsin-related photopigment that may function as 
photopigment or photoisomerase [56], indicating that 
these cells (or adjacent cells) are light-sensitive. Com- 
plementing this, we observed that the previously char- 
acterized c-opsinl [35], an ortholog of rhodopsin and other 
C'Opsins, was expressed in cells around the apical organ 
(Figure 6L). As c-opsinl expression is difficult to score in 
48 hpf larvae (the stage for which the PrImR resource is 
available), we were not able to further characterize the c- 
opsinl-\- cells. One of the otp + cells was demarcated by the 
expression of Imxlab (a terminal selector gene for seroto- 
ninergic neurons [57]) and by position correlates to the 
serotonergic interneuron (Additional file 1: Figure S5A; 
Figure 61). The ampullary apical tuft cells (Additional 
file 1: Figure S4F) specifically expressed hoxl (Figure 2L) 
and a trpV channel previously implicated to serve mechan- 
osensory roles in other protostomes (Figure 2P) [58], 
indicative of multiple sensory modalities. Finally, and 
consistent with their specialized morphology, the two 
putative mechanosensory cells (Additional file 1: Figure 
S5H) expressed miR-183 (Additional file 1: Figure SI), 
a conserved microRNA that demarcates chemo- and 
mechanosensory cells across bilaterians [22,59,60]. Fur- 
thermore, they were the only cells to express otx (compare 
Figure 2F). 

Hierarchical clustering reveals distinct groups of apical 
organ cells 

In line with our observation that at 24 hpf the apical 
organ region including the tuft cells was devoid of six3 
expression, the tuft cells themselves and the directly 
ventrally adjacent sensory-neurosecretory and serotonin- 
ergic cells were 5/;v3-negative at 48 hpf. By contrast, the 
crescent and mechanoreceptor cells expressed six3 and 
thus, by molecular identity, appear to represent differenti- 
ated cell types of the surrounding apical plate. Hierarchical 
clustering based on the 48 hpf PrImR data (Figure 6M) 
supports this distinction: one cluster comprised the flask- 
shaped sensory-neurosecretory cells and the adjacent sero- 
toninergic cell, which, in addition to the absence of sixS 
and presence of otp, miR-7 and phc2 (see above), expressed 
carboxypeptidase-E {cpe) and other genes not directly 
linked to neurosecretion, such as smad2/3. This cluster 
corresponded to the otp + cells devoid of six3 expression at 
24 hpf (see above). Another well-supported cluster 
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comprised the crescent cells and the putative mechanosen- 
sory cells, which, in addition to six3, expressed the micro- 
RNA miR-29 (Figure 6H), the ctbll and bZIRTF genes 
encoding transcription factors of unknown function, and 
teI<tm-2, a structural component of microtubules [61]. 
Note that in our clustering, absence and presence of six3 
and foxq2 strictly correlated, indicative of co-regulation 
[30] (Figure 6M). 

The ampullary tuft cells devoid of six3 (and foxq2) exhib- 
ited a more separate identity with equal distance to either 
cluster. With the six + (and foxq2-\-) mechanosensory and 
crescent cells, they shared expression of foxj, which plays a 
role in the formation of long motile cilia [62], and of three 
microRNAs of the ciliary' group that have previously 
been found to demarcate locomotory ciliary bands across 
Bilateria [22]. They shared the expression of the transcrip- 
tion chxlO with the peptidergic cells (Figure 6M). 

Discussion 

Evolutionary conservation of apical patterning 

We detected strong similarities in molecular topography 
around the apical organ of ciliated larvae in protostomes, 
deuterostomes and cnidarians that we interpret as com- 
mon heritage. It is highly unlikely that a specific and 



partially nested pattern involving at least nine transcrip- 
tion factors {rx, sixS, foxq2, foxj, otx, otp, irx, nkxS, hox) 
evolves convergently twice. The core of this pattern that is 
shared between all neuralians appears to be the apical co- 
expression of two factors, sixS and foxq2, in sea urchin 
[26,27], hemichordate [25], annelid (our data), brachiopod 
[63], sea anemone [28,30] and hydrozoan [64]. In sea ur- 
chin and spiralian larvae, these factors demarcate the 
apical plate and we propose to expand this term to the 
six3-\-,foxq2-\- region of all neuralian primary larvae, which 
should be universally referred to as the apical plate (purple 
in Figure 7). More peripherally, the apical plate co- 
expresses rx, as demonstrated here for the annelid but as 
also found in the sea anemone [28] and hemichordate [25] 
(dark purple in Figure 7). (In sea urchin, rx appears to be 
expressed in the entire apical plate.) In all neuralians in- 
vestigated, the specification of the apical plate is sensitive 
to Wnt signaling, in that pharmacological activation of 
Wnt signaling abolishes the expression of six3, foxq2 
and of other apical markers ([27,28,43,48] and our data). 
The gene regulatory network establishing six3 and foxq2 
expression is only beginning to be elucidated and appears 
to commonly involve activation of foxq2 by sixS [30]. In 
bilaterians, the apical plate also appears to peripherally co- 
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Figure 7 Larval molecular territories and cell types show global conservation of ancient larval patterning. For animals in wliicli molecular 
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express fezf (in annelid, sea urchin and in Saccoglossus) 
and ventrally nkll and hnfS/onecut (in annelid and in sea 
urchin; blue in Figure 7). 

Notably, although the six3-\-, foxq2-\-, rx-\-, fezf-^-, nk2.1, 
onecut + apical plate, as defined here, is a larval charac- 
ter, the same genes are also co-expressed in bilaterian 
phyla that have lost the primary larva, such as insects 
and vertebrates, where they specify anterior brain re- 
gions [34,91,92]. This is consistent with previous obser- 
vations that, after Platynereis metamorphosis, apical 
plate markers remain expressed and demarcate the de- 
veloping cerebral ganglia of the adult nervous system 
[12]. Therefore, while the stereotypical, partially overlap- 
ping co-expression of the above-mentioned genes can be 
used to topographically align' larval body plans, it is not 
a unique feature of primary ciliated larvae. 

A universal molecular definition of the apical organ 
region 

The apical organ develops in the center of the apical 
plate. In annelid, sea urchin [37] and sea anemone [28,30], 
this apical organ region is specifically excluded from the 
region of sixS expression (in annelid (Figure 3 A) and cni- 
darian (Figure 3C) [30]). Instead, it expresses the tran- 
scription factors foxj (in sea anemone [30] and sea urchin 
[38]), nkxS (in annelid (Figure 2Q) and sea anemone [28]), 
irx (in annelid (Figure 2M) and sea anemone [28]) and a 
hox paralog (in annelid (Figure 2L), mollusk [93] and sea 
anemone [28]). Note that all of these factors are also 
expressed elsewhere; for example, in Platynereis^ foxj is 
expressed in other sixS^- ciliated apical plate cells and, 
more generally, in ciliary bands; nkx3 is expressed in other 
apical plate cells (our data) and in segmented mesoderm 
[94] and hoxl is expressed in the second larval segment 
[95] - nevertheless, the recurrent appearance of these fac- 
tors in the apical organ region across neuralians appears 
highly significant and we propose that it reflects the evolu- 
tionary conservation of apical organ cell types (see below). 
In addition, apical organ formation appears to similarly 
depend on local FGF signaling, as shown for sea anemone 
[14] and as suggested by localization of the receptor to 
the apical plate, in sea urchin [37] and annelid (Figure 2N). 
Also, the common localization of Tgf|3 signaling inhibitor 
noggin to the apical organ, as seen in sea anemone [40], 
sea urchin and in Platynereis ^ is indicative of conserved 
signaling events. Together, these transcription factors and 
signaling molecules provide a highly characteristic mo- 
lecular signature for the apical organ region (and the cell 
types therein, see below). Given that the establishment of 
the six3']xo\Q spatially correlates with, and has been ftinc- 
tionally linked to, the formation of the apical tuft (see 
below), we consider this signature a characteristic feature 
of primary larvae and its apparent evolutionary conserva- 
tion lends strong support to the notion that these larvae 



represent an ancient feature of the metazoan life cycle 
(Figure lA and see above). It will be interesting to deter- 
mine to what extent this signature or parts of this signa- 
ture are present in groups that have lost primary larvae; 
given the spotted appearance of irx^ nkxS, hox and foxj 
expression in the apical organ region it is likely that 
any conservation of this signature at adult stages would 
relate to the persistence of apical' cell types throughout 
the life cycle. 

It should be stressed that the hox genes expressed in 
the apical organ in cnidarian, mollusk and annelid repre- 
sent different paralogs of the Hox cluster. Anthoxl, which 
demarcates the apical organ in sea anemone, is a posterior' 
hox gene [96]; LoxS, Lox4 and Lox2, which are expressed 
in different cells of the apical organ of the snail trocho- 
phore [93], belong to the middle part of the cluster; and 
hoxl, that we find expressed in the Platynereis apical 
organ, is an anterior' hox gene. The utilization of different 
hox paralogs in the apical organ can be explained if one as- 
sumes that hox expression in apical organ cells is older 
than the hox cluster itself and already occurred at the times 
when a single m-hox/parahox gene existed [97]. Then, 
concomitant with subsequent duplication events, the ex- 
pression of different hox paralogs would have been lost in 
divergent evolutionary lines. Alternately, the hox paralogs 
utilized in the apical organ may have been exchanged over 
time. Additional studies examining the downstream targets 
of hox expression will be useftil in distinguishing between 
these scenarios. 

The apical tuft - ancient nucleation center of apical 
organs 

If the neuralian ancestors specifically expressed the mo- 
lecular signature defined above in selected cells of the 
apical organ region, what cell types did the partaking 
genes specify and what can we learn about the ancestral 
composition and function of the neuralian apical organ? 
Except for Platynereis, cellular resolution molecular fin- 
gerprints are not yet available in marine primary larvae, 
limiting the extent of interphyletic comparison of specific 
cell types. However, limited molecular evidence is avail- 
able and adds to a rich body of morphological evidence. 
The first and most obvious candidates for conserved ap- 
ical organ cell types are the apical tuft cells that are wide- 
spread in Neuralia [2,11,85,98]. In Platynereis, the apical 
tuft cell exhibits a basket-shaped morphology with charac- 
teristic intracellular support structures resembling those 
of ampullary tuft cells in moUusks [11] and stands out 
with the specific expression of hoxl. In the snail Gibbula, 
expression of Lox4 is confined to 'two cilia-bearing cells of 
the apical organ' [93]; it is as yet unclear how these relate 
to the ampullary tuft cells. In sea anemone, injection of 
anthoxl morpholino results in planulae lacking the apical 
tuft [30], corroborating a link between hox expression and 
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tuft formation. Given that anthoxl knockdown also abol- 
ishes the 'hole' in sixS expression, one might speculate 
that the formation of an apical tuft involving restricted 
hox expression, FGF signaling and apical down-regulation 
of six3 represents the conserved core of apical organs in 
Neuralia (see above). This notion would imply, however, a 
transition from monociliated to multiciliated tuft cells: 
while the Platynereis apical tuft cell bears multiple long 
motile cilia, those of some other spiralian and of deutero- 
stome larvae are monociliated [74]; likewise, the cni- 
darian apical tuft is composed of many cilia that each 
emerge from single monociliated cells [5]. These long 
cilia perform a mechanosensory function during sweeping 
behavior of the benthos [99,100]. In the absence of cell 
type-specific expression data from other neruralian larvae, 
the evolutionary relationship of the various tuft cells re- 
mains speculative; future studies will further clarify their 
evolutionary relationships within neuralians. 

The apical plexus - a conserved sensory-neurosecretory 
release site for the control of body physiology and 
metamorphosis 

Besides the apical tuft, the apical plexus directly under- 
neath the apical tuft may belong to the ground pattern 
of apical organ structures in Neuralia, as it is commonly 
innervated by flask-shaped, neuropeptidergic sensory- 
neurosecretory cells from the apical plate or apical organ 
proper - in mollusks [101,102], sea urchins [78], hemichor- 
dates [74] and, likewise, cnidarians [85,103-107]. Among 
the neuropeptides secreted by these cells into the apical 
plexus may be two ancient types of short amidated pep- 
tides, Wamides and R(F/Y) amides, recently proposed to 
predate neuralians [108]. In cnidarians, RFamidergic sen- 
sory neurons from anterior' larval regions project into the 
apical plexus [85,105,106,109]. Likewise, RFamidergic cells 
populate the apical plate and organ in annelids (as shown 
here and in [12,110]), mollusks [101,102], polyclad flat- 
worms [6] and phoronids [111]. Wamides have been im- 
plicated in the control of larval settlement in cnidarians 
[112,113] and, recently, in annelids [15] and are likewise 
released from the apical plexus. In the late planula larva of 
the hydrozoan Hydr actinia, GLWamide-positive sensory 
cells that innervate the apical plexus populate a belt 
around the apical pole [114], which, by position, forms 
part of the apical plate; MIP, a related Wamide neuropep- 
tide that induces larval settlement in annelids, is secreted 
into the apical plexus by apical organ cells [15]. 

In Platynereis, the transcription factor otp demarcates 
the neuropeptidergic sensory-neurosecretory cells in ap- 
ical organ (Figure 2G,G',G" and [12]). Since otp + cells also 
exist in the apical organ of the snail Patella [84] and are 
overlapping with sixS expression in the apical plate of the 
late Nematostella planula [28], we hypothesize that otp + 
sensory-neurosecretory cells projecting to the apical 



plexus formed part of the apical plate in neuralian an- 
cestors and that subsets of these cells became part of the 
apical organ in the bilaterian lineage. Another gene de- 
marcating the neuropeptidergic sensory-neurosecretory 
cells, among other cells, is nkx3. The apparent deep con- 
servation of nkx3 expression in apical organs (see above) 
suggests an ancient role of nkx3 in the specification of 
these apical neurons; to find out, a thorough investigation 
of nkx3-\- cells in other apical organs, for example in the 
cnidarian planula, will be rewarding. 

Besides the neuropeptidergic cells, apical serotoninergic 
cells represent a third type of neurosecretory cell project- 
ing into the apical plexus. They exist in mollusks [11], 
echinoderms [78], enteropneusts [98] and other bilaterian 
phyla [115]. Serotonergic cells are also enriched in apical 
body regions of the hydrozoan cnidarian Phialidium gre- 
garium, where the release of serotonin has been reported 
to trigger larval settlement [116]. The molecular identity 
of serotoninergic apical cells is beginning to be elucidated: 
the 5HT cell included in our PrImR analysis appears to 
specifically express Imxlah, a LIM homeodomain factor 
implicated in serotoninergic specification in nematodes 
[57], and nk2.1, demarcating serotoninergic cells in sea ur- 
chin [36] and, possibly, in the Ptychodera tornaria [75]. 
Likewise, fezf has been proposed to correlate with seroto- 
ninergic fate in the sea urchin [32]. Further genetic studies 
on the serotonergic system in various marine larvae will 
be needed to resolve this issue. 

Ambient light detection 

Finally, our study provides strong evidence for photo- 
sensitivity being an ancient feature of apical organs. Re- 
markably, the opsins identified in cells in and around 
the apical organ in Platynereis (Figure 6K,L), and also in 
Terebratalia transversa [117] and Nematostella vectensis 
[86], all fall within the peropsin/ciliary opsin families. 
These observations indicate that apical organs evolved 
as multimodal sensory structures, of which photosensi- 
tivity formed a key component. 

'Minimally indirect development' links apical patterning 
of larval and adult stages 

The continuous deployment of the apical patterning sys- 
tem at larval and adult stages and the persistence of some 
apical plate and organ cell types into post-metamorphic 
stages would suggest that a gradual type of metamor- 
phosis (where the bulk of tissues persist, see below) is 
more ancient than the catastrophic' mode of metamor- 
phosis dubbed maximal indirect development [19]. We 
refer to such a biphasic life cycle, with gradual and limited 
metamorphosis in which larval neural structures are in- 
corporated into the adult nervous system, as minimally 
indirect development.' Illustrating this, the Platynereis ap- 
ical organ tuft cell appears to form a nucleation center' 
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around which the brain is organized (Figure 61), and the 
larval axons pioneer the tracts and nerves of the later ner- 
vous system [31,54]. It is possible that the eumetazoan 
common ancestor showed minimally indirect develop- 
ment with a larval stage resembling the primary ciliary lar- 
vae of modern marine bilaterians and cnidarians. 

Conclusions 

We have investigated regionalization of the larval epi- 
sphere, the effects of ectopic activation of Wnt signaling 
on apical patterning, and the molecular fingerprint of ap- 
ical cell types in the marine annelid Platynereis dumerilii. 
Comparing our findings to those in other marine larvae, 
we present a core set of characteristics common to 
primary ciliated larvae in bilaterians and cnidarians. 
All larvae develop an apical plate that we define by a 
combination of transcription factors most prominently 
involving six3 and foxq2. Expression of these factors and 
formation of the apical plate is sensitive to Wnt signaling 
activity. Finally, a conspicuous apical tuft forms within a 
central six3-ivQQ territory within the apical plate. These 
similarities are most parsimoniously explained by com- 
mon origin. We accordingly propose that the last com- 
mon ancestor of bilaterians and cnidarians developed 
via primary larvae that possessed an apical tuft as part of a 
simple apical organ. A basal plexus may have formed dir- 
ectly underneath the apical organ, which was innervated by 
sensory-neurosecretory apical plate cells. We hypothesize 
that an ancient function of the apical organ was the con- 
trol of metamorphosis and opsin-based ambient light 
perception. Various types of additional apical plate cells 
would then have subsequently been recruited to form part 
of the apical organ in the divergent bilaterian lineages. 
Our findings support an ancient and common origin of 
primary ciliated larvae. 

Methods 

Isolation of Platynereis genes and sequence analysis 

Partial sequences for pdu-fgfr, pdu-foxq2, pdu-irx, pdu- 
foxjl, pdu-hnf6, pdu-wnt4, pdu-frizzled5/8, pdu-sfrpl/5 
and pdu-peropsin were assembled from Platynereis tran- 
scriptome and genome resources, amplified with spe- 
cific primers, and extended with rapid amplification 
of cDNA ends PGR. Pdu-CTBLl, Pdu-bZIP-TF, Pdu- 
tektin-2, Pdu-spondin, Pdu-gbrl, Pdu-cpe, Pdu-smad2/3 
and Pdu-klf2/4 were characterized during a whole mount 
in situ hybridization screen from expressed sequence tag 
clones (RT and DA, unpublished). The GenBank ac- 
cession numbers for peropsin, foxj, fezf, onecut, fgfR, 
noggin, foxq2, frizzled4, frizzled5/8, frizzled9/10 and 
sfrpl/5 are [GenBank:KF844232] to [GenBank:KF844242] 
respectively. Accession numbers for ctbll, spondin, gbrl, 
cpe, ces2 and klf are [GenBank:KF835846] to [GenBank: 
KF835851] respectively. 



In situ hybridizations and immunostainings 

For in situ hybridizations at early stages, Platynereis 
larvae were fixed in 4% PFA, 0.1 M MOPS, 2 mM EGTA, 
1 (iM MgS04 and 0.1% Tween-20, for 4 to 6 hours at 4°C, 
then rinsed in PTW and ice cold methanol, followed by 
storage in methanol at -20C. 

In situ hybridization for early Platynereis larvae were 
performed as in [118,119] with the following modifications. 
Embryos were digested with 0.1 mg/mL proteinaseK for 
30 seconds (Merck 7066304, Darmstadt, Germany). Fol- 
lowing hybridization, 0.5 x SSC washes were replaced by 
0.15 X SSC washes for 15 rather than 30 minutes. In situ 
hybridizations were performed for 48 hpf Platynereis and 
microRNA as previously published [22,118]. The axonal 
scaffold was counterstained with an antibody against 
tyrosinated or acetylated-tubulin (1:500, Sigma T9028 
and T6793). Immunostainings were performed as de- 
scribed [12] with the following primary antibodies: mouse 
anti-acetylated tubulin (1:500, Sigma T6793), rabbit anti- 
serotonin (1:500, ImmunoStar 20080) and rabbit anti- 
FMRF (1:200, Phoenix Pharmaceuticals, H-047-29). 

For the staining of mitotic cells, Platynereis larvae 
were incubated in 10 (iM EdU (Click-iT EdU Imaging 
Kit, Invitrogen, C10340) from 22 to 24 hpf; EdU incorp- 
oration was detected after the incubation with secondary 
antibodies, following manufacturer instructions. 

Alsterpaullone and azakenpaullone treatments 

Platynereis larvae were incubated from 12 to 24 hpf in 
five different concentrations of azakenpaullone (0.1 to 
10 \M) in 0.5% dimethyl sulfoxide (DMSO) in fiftered 
seawater. The number of embryos displaying wild-type, 
reduced or expanded expression patterns for episphere 
molecular markers were assessed from two different bio- 
logical replicates (numbers of embryos displaying each 
pattern are indicated in Figure 5). As most genes assayed 
were significantly affected at 1 and 5 (iM concentrations, 
we conducted washout experiments at these concentra- 
tions. Larvae were exposed to azakenpaullone from 12 
to 24 hpf and subsequently washed out of the pharmaco- 
logical treatment into 0.5% DMSO and incubated in 
0.5% DMSO from 24 to 30 hpf. Controls were main- 
tained in 0.5% DMSO from 12 to 24 hpf and washed 
into a new 0.5% DMSO treatment from 24 to 30 hpf. 
Embryos were examined and assessed in relation to 
control-treated embryos. Effects of alsterpaullone were also 
assayed, compared to those of azakenpaullone, and found 
to be similar to azakenpaullone treatments (Additional 
file 1: Figure S2). 

Microscopy and image processing 

We used reflection microscopy [120] to acquire confocal 
images of in situ hybridization stainings. Fluorescent 
signals in Figure 2Q', Figure 4C and Figure 4C' were 
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obtained using the fluorescent signal emitted from NBT/ 
BCIP precipitate, an alternate method to reflection micros- 
copy. Confocal stacks were tal<en on a Leica TCS SPE with 
a 40x oil immersion objective. Images were processed with 
ImageJ, utilizing either brightest point or average intensity 
settings to generate projections. Subsequently, images were 
cropped and processed in Photoshop, Adobe, San Jose, 
California, USA; brightness and contrast were adjusted 
equally across the entire image. 

For PrImR, average expression patterns at 48 hpf were ob- 
tained after image registration of in situ confocal scans on a 
common reference axonal scaffold, as described in [47]. 

Molecular fingerprint analysis 

Morphologically distinct apical organ ceU types were 
identified by analyzing immunostained larvae at early 
stages (24 and 30 hpf). These ceUs were then located in 
48 hpf larvae for the gene expression analysis with PrImR. 
The co-localization between two average gene expression 
patterns was inspected and visualized using the 'Colocali- 
zation highlighter' plugin in ImageJ. 

Whenever a PrImR average expression pattern was not 
available for the gene and/or the stage of interest, speci- 
mens stained with the gene of interest and tyrosinated 
tubulin were inspected under fluorescence microscopy. 

Hierarchical clustering of molecular fingerprints was 
carried in R (the R Project for Statistical Computing) 
from the dataset in Additional file 1: Figure S5, using 
Pearson correlation and average linkage. 

Phylogenetic analyses for gene orthologies 

Platynereis dumerilii gene coding sequences used in this 
study were isolated as described above. Sequence data from 
the lophotrochozoans Lottia gigantea and Capitella teleta 
and the cnidarian Nematostella vectensis were identified on 
their respective JGI genome portal webservers. Additional 
sequences used in the analyses were downloaded from 
Treefam [121]. Multiple alignments of predicted proteins 
were generated with MUSCLE using the default settings 
[122] and were subsequently inspected and corrected by 
eye. Full alignments were trimmed using G-blocks [123] 
and were run through ProtTest using the default settings 
to determine the optimum evolutionary model for phylo- 
genetic analyses [124]. Neighbor joining trees were con- 
structed using MEGA [125] and maximum likelihood 
analyses were conducted using PhyML [126] with the 
amino acid substitution models specified from ProtTest. 
Phylogenetic trees are available in Additional file 2. 

Additional files 



Additional file 1: Figure SI. miR183 expression in a 48 hpf larva. 

Figure S2. Alsterpaullone treated embryos. Figure S3. Azakenpaullone 
treated embryos following washout. Figure S4. Percentages of affected 



embryos in azakenpaullone treatments. Figure S5. PrImR analysis of the 
expression of transcription factors, miRNAs and differentiation markers in 
defined apical organ cell types. 

Additional file 2: Phylogenetic analysis for the assignment of gene 
orthology. 
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